A BOUNDARY POINT LEMMA FOR BLACK-SCHOLES 
TYPE OPERATORS 
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| Abstract. We prove a sharp version of the Hopf boundary point lemma 

f*^ . for Black-Scholes type equations. We also investigate the existence and 

' the regularity of the spatial derivative of the solutions at the spatial 

^ ' boundary. 

^ : 
m . 

^ \ 1. Introduction 

,— ,- In this paper we study a parabolic initial-boundary value problem mo- 

Ph \ tivated by applications involving modeling of non-negative quantities de- 

•<r^ ■ scribed by diffusion processes that are absorbed at 0. To describe the prob- 
lem in a financial setting, consider a market consisting of a risk-free asset 

d ■ dB = r(t)Bdt 

a- 

where r is a deterministic function, and n risky assets with non-negative 
prices Xi, i = l,...,n, modeled by diffusion processes. According to stan- 
dard arbitrage theory, to price options maturing at some future time T 
■ written on such risky assets, we need not specify the expected rate of re- 

(T) , turn of the assets, compare for example Theorem 7.8 and Proposition 7.9 in 

CN ' [2] and Proposition 2.2.3 in ^Uj. Instead one specifies the price dynamics 

under a so-called risk-neutral measure, under which the expected rate of 
^ | return of all assets equals the interest rate r(t). For the sake of notational 

convenience, however, replacing Xi(t) by Xi(t) exp{J^ r(s) ds}, we obtain 
processes with drift 0. Financially this corresponds to quoting the assets in 
terms of bonds maturing at time T. By abuse of notation, we denote also 
these new processes by X{. Thus we model the price dynamics of the ith 
asset by 

dXi = <*ij (X i} t) dWj Xi(t) = Xi . 

Here Wj, j = 1, n, are independent standard Brownian motions, and a = 
i a ij)?j=i 1S f° r eacn x an d t an nxn-matrix. We assume that the components 
of X = (Xi, ...,X n ) are absorbed at 0, i.e. if some component at some point 
reaches 0, then it remains forever. Therefore we let a;y(0, t) = for all i 
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and j. We also assume that the rank of a is equal to the number of non-zero 
spatial coordinates (this assumption is connected to the completeness of the 
model; more precisely, under this assumption every option has a unique 
arbitrage-free price, compare Corollary 12.2.6 in 13 or Theorem 1.6.6 in 
jl()j). Further assumptions on a are specified below. Given a continuous 
contract function g : [0, oo) n — > R of at most polynomial growth, the value 
U{x,t) at time t of an option paying g(X(T)) at time T is 

U{x,t) :=E x>t g(X{T)), 

where the indices indicate that X(t) = x. Alternatively, the function 
u(x,t) := U(x,T — t) solves the Black-Scholes equation 

m j Cu = for (x, t) G (0, oo) n x (0, oo) 

( j \ u(x,0) = g(x), 

where 

n 

Cu = ^2 aijU XiXj - Uf, 
i,i=l 

ciij = a>ij(x,t) are the elements of the matrix 

a{x, t) := a{x, T — t)a*(x, T — t)/2 

and 

du d 2 u 
ut := ^, u XiXj := 

In order to have uniqueness of solutions to the problem in the class of 
functions with at most polynomial growth, one needs (in general) to im- 
pose boundary conditions at the faces {xi = 0}, i = l,...,n. We assume 
that these boundary conditions are defined inductively by solving the par- 
tial differential equation in lower dimensional faces (the assumptions on a 
guarantee that the operator is parabolic in these faces). One thus starts 
with solving an ordinary differential equation (ut = 0) along the t-axis, then 
one solves n parabolic equations in the faces spanned by t and one of the 
variables xi,...,x n , and so on. According to Theorems 4.1 and 5.5 in 
this procedure results in a unique classical solution of polynomial growth to 
equation 

In this paper we are concerned with the behavior of the first order deriva- 
tives u Xi at the boundary. We should note that the derivatives u Xi are of 
special importance in financial applications. They are the so-called "deltas" 
of the option, and they represent the number of stocks Xi a hedger should 
have in a hedging portfolio, compare for example Theorem 8.5 in [2j. 

The following theorem, the Hopf boundary point lemma (adapted to our 
current setting), is well-known in the theory of parabolic equations, compare 
for example 6 (Theorem 2.14, p. 49) or (Lemma 2.6, p. 10). 

Theorem 1.1. (Hopf boundary point lemma.) Assume that the differ- 
ential operator is uniformly parabolic, i.e. that there exists 7 > such that 
£a(x, t)£* > 7|£| 2 for all £ E R n , x and t. Let a point P' = (0, x' 2 , ...,x' n ,t ) 
with x'2, x' n , to > be given. Assume that in a neighborhood of P' we have 
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u(P) > u(P') for all points P with x\ > and t < to- Then u Xl (P') > in 
the sense that 

liminf-(u(P' + eei) - u(P'j) > 0, 



e— o e 

where e\ is the unit vector in the x\-direction. 

The following well-known example shows that the above theorem is not 
valid without the assumption of uniform parabolicity. 

Example (The Black-Scholes price of a call option.) Let n = 1, and 

let a(x, t) = ax and g(x) = (x — K) + where a and K are positive constants 
(here and in the sequel we drop the subscripts of cm, an and x\ if n = 1). 
Then 

«(,,«) = - k<s,^" ( J_' h ' 2 ) 

where 



1 f z 

$(z) = -= exp{— u 2 /2} du. 

V27T J_oo 



It is now straightforward to check that the delta of the call option is given 
by 

'\n(x/K)+a 2 t/2^ 



u x {x,t) = $(- 



vVt 
so u x (0, t) = 0. 

The outline of the present article is as follows. In Section [21 we provide a 
Hopf lemma for equations that are not uniformly parabolic, compare Theo- 
rems and E3I To prove these results we need to assume that an > Cx^ 
for some < (3 < 2. In view of the above example (in which an = a 2 xf/2), 
our result can be regarded as a sharp version of the Hopf lemma. Along the 
lines of the above example we also show that if n = 1 and a < Cx 2 , then 
the result of the Hopf lemma always fails, compare Theorem 12.31 

In Section El we perform further investigations of the spatial derivative at 
the boundary. We use preservation of convexity in dimension one to show 
that u x (0,t) always exists finitely for a wide class of initial conditions, and 
we provide an example in which u x (0,t) is discontinuous. 

2. Hopf boundary point lemma 

In this section we prove a version of the Hopf boundary point lemma for 
non-uniformly parabolic operators. 

Hypothesis 2.1. The nxn-matrix a = (ctij(x,t))™ ._ 1 is defined on [0, oo)" x 
(— oo,T], and 

(i) is continuous in the time variable, and a is also Lipschitz continuous 
in the spatial variables on every compact subset o/(0, oo)™ x (— oo, T]; 
we also assume the standard growth condition \a(x,t)\ < C(l + |x|) 
for some constant C; 

(ii) for each pair (i,j), otij is a function merely of Xi and t (this is 
automatic if n = 1), and for Xi = we have ctij = 0; 

(iii) the rank of a is equal to the number of non-zero spatial coordinates. 
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Ifn = 1 the local Lipschitz condition may be replaced with a local H6lder(l/2) 
condition. 

Note that we allow discontinuities of the coefficients at the spatial bound- 
ary. Also note that condition (ii) does of course not rule out the possibility 
of dependence between the different assets. It merely means that the in- 
stantaneous covariance a™ between two assets Xi and Xj only depends on 
time and the present values of Xi and Xj. The condition (hi) is connected 
to the issue of completeness of the model, compare the discussion in the 
introduction. 

Theorem 2.2. Assume Hypothesis \2.1\ that u satisfies the Black-Scholes 
equation £JJ), and that an(x\,t) > Cxf for some constants C > and 
(3 € [0,2). Let a point P' = (0, x' 2 , x' n , to) with to > and x\ > for 
i = 2, ...,n be given. Assume that in some neighborhood of P' we have that 

u(xi,x 2 , -,x n ,t) > u(0,x 2 ,...,x n ,t) 

for all points with x\ > and t < to. Then u Xl (P') > (in the same sense 
as described in Theorem 

Proof. Introduce the function 

n 

v(x, t) := xi + x\ +e -\t- t \ N ~Y,\xi- x'^ 

1=2 

for some constants e > and 2V > 1, both to be chosen later. Let 

D:={P : v(P) > 0,0 < x x < rj,t < t } 

for some small constant 7] > 0. Then there exists a positive constant C\ 
(depending on r\) such that 

-C lX \ /N <t-t <0 

and 

-C lX l /i2N) < -\xi-xH < 0, 
i = 2, ...n, for points in D. Consequently, 

n n 

£v = aijV XiXj - Vt=/ J auVxiXi ~ v t 

i,j=l i=l 

n 

= an(l + e)ex\- 1 - N\t - tof" 1 - (AN 2 - 2N) £ a ti \xi - x\\ 2N - 2 

i=2 

> C (l + e)ex^- 1 - NC 2 x[ N - 1),N - (4iV 2 - 2N)C 2 xf- 1)/N 
in D for some constant C 2 satisfying 

C 2 > max {Cf-\cf- 2 (n - 1) max {a u (P')}}. 

2<i<n 

Thus, choosing e small and N large so that (3 + e — 1 < (N — 1) /N < 1, it is 
clear that Cv > in D (at least if rj is small enough; note that C\, and thus 
also C 2 , can be held fixed when decreasing rf). The parabolic boundary d p D 
of D can be written as d p D = S\ U S 2 where 

S 1 = {P :v{P) =0,t<t ,x 1 <7 ] } 
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and 

S 2 = {P :v(P)>0,t<t ,x 1 =r ] }. 

Since, by assumption, 

u{x x ,x 2 ,...,x n ,t) > u(0,x 2 ,...,x n ,t) for (x\, x 2 , -, x n , t) G D\{P'} 

there exists 5 > such that u(x±, x 2 , x n , t) — u(0, x 2 , x n , t) > 5 for 
(x±,x 2 , ...,x n ,t) G S 2 . Since v is bounded by 1 in D (at least if r] is small 
enough) it follows that 

(2) u(xi,x 2 , ...,x n ,t) - u(0, x 2 , ...,x n ,t) - 5v(xi,x 2 , ...,x n ,t) > on d p D. 
Moreover, applying the differential operator 

d 2 d 



i,j=l 



'j 



dxidxj dt 



in D to the function u(0, x 2 , x n , t) we get 

(3) C\u(0,x 2 ,...,x n ,t)j = y~] aij(xi,Xj,t) u(0, x 2 , x n ,t) 

i,j=2 •'' 

-—u(0,x 2 , ...,x n ,t) = 

since n(0, x 2 , ...,x n ,t) satisfies the (n — l)-dimensional Black-Scholes equa- 
tion in the face x\ = 0. Thus we have 

(4) c(u(x!,x 2l ...,x n ,t)-u(0,x 2 ,...,x n ,t)-5v(xi,x 2 ,...,x n ,tj^ < in D. 

Applying the weak maximum principle to the inequalities (J2J) and Q we 
find that 

u(x 1 ,x 2 , ...,x n ,t) - u(0,x 2 , ... ,x n ,t) - 8v(xi,x 2 ,...,x n ,t) >0'mD. 
Using v(P') = we get 

u Xl (P')>5v Xl (P')=5>0 
which finishes the proof. □ 

Remark Note that the assumption that Ojj is a function merely of Xi and 
t (condition (ii) in Hypothesis 12.11) is essential. Indeed, if a instead would 
depend on the whole vector x and t, then the equality © would not be true 
in general. 

Remark Theorem 12.21 also holds if one includes lower order terms in the 
differential operator C. In addition to the assumptions on the coefficients 
specified in Theorem 12.41 below, one also needs to assume that c and 6j, 
i = 2,...,n are independent of x\, where c and hi are as defined in that 
theorem. 

The theorem also generalizes to any non-tangential direction with non- 
positive time-component. In that case, N has to be chosen strictly larger 
than 1. 

To illustrate the result of Theorem 12.21 we give two examples. 
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Example (The Margrabe exchange option.) Let n = 2, assume that 

dX i = u i xf i/2 dWi, 

where U{ > and < < 2 for i = 1, 2, and let 5(^1, a^) = (^2 — 
From Theorem 12.21 it follows that u X2 (xi,0,t) > for x 2 > and t > 0. 
To investigate the derivative at the boundary xi = 0, note that adding 
the affine function x\ to the contract function g gives (x2 — x\) + + x\ = 
maxjxi, X2}- By Theorem 12.21 the value u of this contract has a spatial 
derivative u xi (0, a; 2 ,i) > 0. It follows that u xi = ]^(u — x\) > —1 for 
points (0, X2,t). 

On the other hand, if j3\ = @2 = 2, i.e. if X\ and X2 are geomet- 
ric Brownian motions, then standard formulas for the value of the ex- 
change option (compare ,12^ ) can be used to show that u X2 (xi,0,t) = 
and u Xl (0,X2,t) = —1. 

Example (The call option in a CEV-model.) Let n = 1, g{x) = 
(x — K) + and 

dX = aXP /2 dW 

for some constants K > 0, a > and (3 E [0, 2]. If j3 = 2, then X is a 
geometric Brownian motion and we know that u x (0,t) = (compare the 
example in the introduction), and if (3 = 1, then one can use the valuation 
formula in 3 (p. 161, equation (36)) to explicitly calculate that u x (0,t) = 
exp{ — ^t} > for t > 0. Theorem 12.21 tells us that also in the remaining 
cases, i.e. for all < 2, we have u x (0,t) > for t > 0. 

In dimension n = 1 we have the following result that for instance covers 
the example given in the introduction. It shows that the condition that 
an > Cx^l for some f3 G (0, 2) in Theorem 12.21 cannot be substantially 
weakened. 

Theorem 2.3. Let n = 1, and assume that a(x, t) < Cx 2 for all x and t for 
some constant C. Also assume that g' exists at x = 0. Then u x (x,t) exists 
at x = and u x (0,t) = g'(0) for all t. 

Proof. By subtracting a constant and a suitable multiple of x we may with- 
out loss of generality assume that g(0) = and </(0) = (note that all 
affine functions w satisfy Cw = 0). Let e > be given. Since g is of at most 
polynomial growth, we can find Ci > and N > 1 such that 

g(x) < ex + C\x N 

for all x. Note that Ci depends on e, whereas ./V can be held fixed if varying 
e. Define the function v by 

v(x, t) := ex + Cix N + C 2 tx N 

for some constant C2 > to be chosen. Then 

v t - a(x, t)v xx > C 2 x N - Cx 2 CiN(N - l)x N ~ 2 - Cx 2 C 2 tN(N - l)x N ~ 2 

= (C 2 - C(d + C 2 t)N(N - l))x N , 
so for times < t < 2C(N i -N) =: *° we have 

v t - a(x,t)v xx > (C 2 /2 - CCiN(N - l))x N > 
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if C2 is chosen large enough. Thus v is a supersolution to the Black-Scholes 
equation satisfying v(x,0) > g(x). It follows from the maximum principle 
that u < v for all times < t < to- Thus u x (0, t) < v x (0,t) = e for such t. 
The same argument applied to — u gives u x (0,t) > —v x (0,t) = — e. Letting 
e->0we arrive at u x (0, t) = for times t < to. Viewing u(x, to) as the initial 
condition the above argument shows that u x (0, t) = also for to < t < 2to 
(note that the same constant N can be used again) and thus by iteration 
also for all t > 0, which finishes the proof. □ 

There is no immediate generalization of Theorem l2.3l to higher dimensions. 
For instance, the price u(x, t) in a geometric Brownian motion model of a 
call option on the sum of two assets (i.e. g(xi,X2) = (x± + xi — K) + ) 
satisfies u Xl {0, X2,t) 7^ g Xl (0, x^)- On the other hand, the example with the 
exchange option shows that the condition on a\\ in Theorem 12.21 is sharp 
also in several dimensions. 

Note that in the version of the Hopf boundary point lemma provided 
above (Theorem l2.2jl it is not assumed that P' is a minimum of the function 
u but rather of the function u(x\, X2, x n , t) — ii(0, X2, x n , t). This is 
appropriate in the case when the boundary conditions are defined inductively 
by solving the equation in lower dimensional faces. Our proof also works in 
the situation where u(xi,X2, ...,x n ,t) has a minimum in P' . This situation 
is not the typical one when the boundary conditions are defined as in the 
introduction, but it can of course occur in other types of degenerate initial- 
boundary value problems. 

Theorem 2.4. Assume that a is continuous and parabolic in (0, oo) n x 
(0, 00), i.e. £a(x,i)£;* > for all £ G W n \ {0} and for all x,t. Assume also 
that u satisfies 

n n 

(5) ^2 aijU XiXj + ^2 biU Xi + cu-u t < 

i,j=l i=l 

in (0, oo) n x (0, 00), where hi = bi(x, t) and c = c(x, t) are continuous func- 
tions. Let a point P' = (0, x' 2 , x' n , to) with x' 2 , x' n , to > be given, and 
assume that in a neighborhood of P' we have 

an > Cxi, 
an < C 

fori = 2,...,n, 

(6) 61 > -Cx^ 1+S , 

\bi\ < Cxt 2+& 

for i = 2, n and 

(7) c > -Cx^- 2+S 

for some constants C > 0, 5 > and (5 G [0,2). Also, in the same neigh- 
borhood, assume that u{P) > u(P') for all points P with x\ > and t < to- 
Then u Xl (P') > 0. 

Proof. The proof follows along the lines of the proof of Theorem 12.21 and is 
therefore omitted. □ 
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As is shown in the next example, the lower bounds on b and c in Theo- 
rem 1^3] are sharp in the sense that the result fails for 5 = 0. 

Example Let n = 1. The function u(x,t) := x 2 /2 is the solution of 



Note that @ is not satisfied and that u x (0,t) = 0. Moreover, the same 
function u also is the solution of 



\ u(x,0) = x 2 /2. 
For this system Q is not satisfied. 

Remark If (3 = in Theorem 12.41 i.e. if the differential operator is uni- 
formly parabolic, then one only needs the inequality © to hold in a parabolic 
frustrum, compare Lemma 2.6, p. 10 in 11 j. Examining the proof of the 
Hopf lemma for non-uniformly operators, it is clear that this also is true for 
(3 G [0, 1), since N in that case can be chosen to be 1. 



In this section we use preservation of convexity for n = 1 to deduce the 
existence of the derivative u x (0,t) and also some regularity properties of 
this function as a function of time. If the operator is uniformly parabolic, 
then u x (0,t) exists and is continuous for t > (even if <7'(0) does not exist). 
To see this, note first that we can assume, without loss of generality, that 
g(0) = 0. Then, extending g to a continuous odd function on R and a 
to an even function in x, standard interior regularity results for parabolic 
PDE:s yield that u x (0, t) is continuous. Note that, indeed, the solution to 
the extended problem agrees with the solution to the original problem for 
non- negative x, since the extended solution is odd in x and thus vanishes at 



We start this section with an example that shows that this is not true in 
general for degenerate operators. 

Example (Power options in a geometric Brownian motion model.) 

Let a(x, t) = x 2 and g(x) = x 1 for some constant 7 > 0. Then it is easy to 
check that 



Thus u x (0, t) does not exist if < 7 < 1. 

Theorem 3.1. Assume that the contract function g is convex and that 
g'(0+) > —00. Then the derivative u x (0,t) exists for all t. Moreover, the 
function 1 1— > u x (0,t) is increasing and upper semi- continuous. 

Proof. Recall, see [H], that convexity is preserved and that the option price 
increases in time to maturity, i.e. x 1— > u(x, t) is convex for each fixed t > 
and 1 1 — > u(x,t) is increasing for each fixed x, see also [S] and Since 
u(0,t) = g(0) it follows that u x (0,t) exists and that u x (0, t) is increasing. 
Moreover, using the continuity in t and the spatial convexity of u, it follows 
that u x (0,t) is upper semi-continuous. □ 





3. The spatial derivative in dimension one 



x = 0. 



u(x, t) = x J exp{(7 2 — 7)^}. 
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Remark Note that u x (0,t) also exists if the contract function can be writ- 
ten as a difference of two convex functions (both with finite derivative in the 
origin). Also note that the example above with g{x) = x 7 where < 7 < 1 
(or rather g{x) = — x 7 ) shows that the assumption in Theorem 13.11 about 
g'(0) > —00 is essential. 

Remark In higher dimensions convexity is in general no longer preserved, 
compare (!]. Therefore the above proof of Theorem 13.11 is not applicable 
to problems with several underlying assets. It remains an open question 
to determine conditions under which the spatial derivatives exist at the 
boundary if n > 2. 

The conclusions of Theorem 13.11 that can be drawn more or less imme- 
diately from the preservation of convexity cannot be improved in the sense 
that u x (0,t) need not be continuous in t. Indeed, we end this article with 
the construction of an example where a is continuous and locally Lipschitz 
in x, the contract function g is convex and in C°°([0,oo)) and yet u x (0,t) 
fails to be lower semi-continuous. 

Example (u x (0,t) need not be continuous as a function of t). The 

example is constructed by patching together two functions v : [0, 00) x 
[0, to) — ► [0, 00) and w : [0, 00) x [t , 00) — > [0, 00) for some given to > 0. 

To define v, let h £ C°°([0,oo)) be convex, non-negative and satisfy 
h(0) = 0, fc'(0) = 0, and h'(y) = 1 and h"(y) = for y > y > 0. Let 

C := lim (yh'(y) - h(y)) = y - h(y ). 
Define v : [0, 00) x [0, to) — ► [0, 00) by 

v(x,t) = xV + (t - t)h( 



It follows that 



to-*' 
v t = a(x,t)v xx 



in (0, 00) x (0, to) where 

xV(t - t) + xh'(^r t ) - (t - t)h(^) 



a(x, t) 



2e t {t -t) + h"{, x 



Note that 7^(0, t) = for < t < to, and for x > we have 

lim v(x, t) = x 2 e*° + x, 

lim v t (x,t) = xV° + C, 

(y,t)->(x,t ) 

lim v x (x,t) = 2xe*° + 1 
(y,t)-*-(x,to) 



and 



lim v xx (x,t) = 2e*°. 



Next define w as the unique solution to 



if x = 

x 2 e*° + x if t = t . 
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Then it is straightforward to check that 

w t (x,t ) = xV° + C, 



and 

for x > 0. It follows that 

u(x, t) 

solves 
where 

a(x, t) : 



w x (x,t ) = 2xe to + 1 
w xx (x,t ) = 2e to 



v(x,t) if < i < t 
w(x, t) if to < t < oo 

u t = a(x,t)u xx 

a(x,t) if < t < t 
if t < t < oo. 



Note that a := V2a is continuous on (0, oo) x (0, oo), locally Lipschitz in 
the x-variable and that a satisfies the growth condition 

a(x,t) < D{l + x) 

for some positive constant D. Recall that option prices with convex con- 
tract functions are increasing in the time to maturity, compare [Hj, i.e. the 
function t i— ► w(x,t) is increasing. Thus w x (0,t) > 1 for t > to- Therefore 
the function u x (0,t) is not continuous. 
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